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ABSTRACT 


This  report  provides  the  theoretical  analysis,  within  the  framework  of  special 
relativity,  for  the  information  content  of  a  radial  velocity  spectrometer.  This  device 
is  an  automated  telescope-photometer-spectrometer  combination.  It  measures  the 
displacement  between  a  Doppler-shifted  set  of  spectral  lines  and  its  at-rest  version. 
Thus,  it  provides  a  measurement  of  “radial  velocity.”  Exactly  what  kinematical 
content  is  revealed  by  such  a  measurement  within  the  solar  system  is  brought  out 
herein.  The  discussion  is  couched  in  terms  of  a  satellite-based  sensor  observing 
another  orbiting  body.  Both  objects  are  presumed  to  orbit  about  the  Earth. 
Sufficient  generality  is  presented  to  treat  any  other  case  also.  In  the  specific  case  of 
interest,  both  the  topocentric  distance  and  the  topocentric  radial  velocity  are 
present  in  the  measured  quantity. 
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RADIAL  VELOCITY  SPECTROMETERS: 
WHAT  THEY  MEASURE 

I.  INTRODUCTION 


Passive  optical  measurements  of  artificial  satellites  and  ballistic  missiles  are  usually  thought 
of  as  providing  angles-only  information;  that  is,  such  observations  can  only  provide  the  direction 
to  (or  position  of)  a  source  of  light.  As  should  be  well  known  from  the  astronomical  literature, 
much  more  can  be  learned  about  the  light  source’s  nature,  center  of  mass  motion(s),  and  dynami¬ 
cal  state.  In  particular,  by  observing  the  shift  of  various  spectral  lines  from  their  at-rest  locations 
one  can  determine  the  radial  velocity  of  the  source.  For  stars  and  galaxies  the  nonrelativistic 
approximation  may  be  written  as 

AX/X  =  vr/c 

wherein  AX  is  the  signed  wavelength  shift  at  wavelength  X  (i.e. ,  the  Doppler  shift),  c  is  the  speed 
of  light  in  vacuo,  and  vr  is  the  signed  component  of  the  projection  of  the  source’s  velocity  vector 
along  the  line  of  sight  to  the  source  (i.e.,  the  radial  velocity).  The  physics  of  the  Doppler  shift 
applies  equally  well  to  luminous  bodies  and  to  those  visible  owing  to  reflected  light.  Thus,  in 
principle,  artificial  satellites  too  may  have  their  Doppler  shifts  measured  by  passive,  optical 
means. 

However,  in  the  case  of  artificial  satellites,  the  dynamical  situation  is  a  little  more  complex 
than  it  is  for  stars  or  galaxies.  Indeed,  the  above  formula  is  an  oversimplification  for  these  celes¬ 
tial  objects  too  as  the  telescope  is  (usually)  on  the  surface  of  a  rotating  Earth,  the  Earth’s  center 
of  mass  is  in  motion  about  the  Earth-Moon  barycenter,  the  Earth-Moon  barycenter  is  in  motion 
about  the  Sun,  and  the  Sun  is  in  motion  about  the  center  of  mass  of  the  solar  system.  The  solar 
system  barycentric  reference  frame  is  the  best  approximation  that  we  have  to  a  realizable  inertial 
frame.  (For  galaxies  one  should  add  yet  other  components,  as  the  solar  system  barycenter  is  in 
motion  relative  to  the  nearby  stars  and  the  center  of  the  Milky  Way.) 

All  the  corrections  necessary  owing  to  the  conflict  between  the  natural  reference  frame  of  the 
observations  and  the  desire  to  be  in  an  inertial  reference  frame  are  very  small  and  simple  to  per¬ 
form  nonrelativistically.  The  typical  precision  in  galaxy  work  is  —  50  km/s,  while  for  stellar  work 
1  km/s  is  considered  good.  Since  the  entire  range  of  satellite  speeds  is  only  —  20  km/s,  a  much 
higher  precision  is  required  in  order  to  be  a  useful  addition  to  the  endeavors  of  initial  orbit  deter¬ 
mination,  rapid  maneuver  detection,  or  to  the  highly  refined  differential  correction  of  orbits.  This 
not  only  puts  severe  demands  on  the  instrumentation,  but  requires  a  more  sophisticated  analysis 
of  exactly  what  is  measured.  That  analysis  is  contained  in  this  Report. 

Before  getting  to  the  analysis,  a  few  words  about  the  potential  instrument  are  in  order.  The 
classical  astronomical  technique  of  radial  velocity  determination  utilizes  a  spectrograph  to  photo¬ 
graph  the  dispersed  light  from  the  star  and  then  to  measure  the  line  shifts  after  developing  the 


photographic  plate.  This  is  a  slow  and  very  time-consuming  process,  especially  as  all  the  light 
from  the  source  is  being  used  to  ascertain  only  one  quantity.  Moreover,  this  single  number  repre¬ 
sents  a  global  aspect  of  the  entire  spectrum.  Thus,  first  Fellgett*  and  then,  especially,  Griffin+ 
proposed  utilizing  a  cross-correlation  technique  at  the  telescope  to  increase  the  efficiency  of  the 
measurement.  Such  automated  devices  are  now  known  as  radial  velocity  spectrometers.  Their 
speed  of  operation  is  such  that  a  precision  of  —  10  m/s  on  a  10m  object  during  an  integration 
time  of  —  10  s  seems  feasible.  Future  reports  shall  dwell  on  the  design  of  such  a  device  keyed  to 
the  solar  spectrum. 


*  P.B.  Fellgett,  Optica  Acta  2,  9  (1953). 
t  R.F.  Griffin,  Astrophys.  J.  148,  465  (1967). 
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II.  COMPLICATIONS 


The  artificial  satellites  or  ballistic  missiles  that  one  might  wish  to  observe  with  a  radial  veloc¬ 
ity  spectrometer  shine  by  reflected  sunlight.  Therefore,  the  shape  of  the  line  spectrum  of  the 
Sun’s  light,  the  temporal  variation  of  these  features,  and  the  nature  of  the  continuum’s  envelope 
are  of  great  importance.  The  Sun  can  be  regarded  as  a  point  source;  thus,  such  problems  are 
greatly  alleviated  because  one  is  effectively  integrating  over  the  whole  —  32'  solar  disk.  This 
washes  out  much  of  the  rapid  temporal  and  small-scale  spatial  variability.  Also,  the  Sun  is  a  very 
well  studied  object.  This  fund  of  knowledge  can  be  utilized  to  minimize  data  reduction  problems. 
In  particular,  one  can  concentrate  on  that  portion  of  the  solar  spectrum  arising  from  the  more 
stable  photosphere  rather  than  from  the  less  stable  chromosphere  or  corona.  A  more  detailed  dis¬ 
cussion  of  these  points  will  appear  in  a  future  Report,  and  herein  I  shall  merely  quote  the  essen¬ 
tial  point;  microturbulence  in  the  solar  atmosphere  occurs  on  the  —  3  m/s  level  on  time  scales 
short  enough  to  limit  one’s  ultimate  precision  to  this.  Note  that  this  is  about  a  part  in  ten  thou¬ 
sand  of  the  quantity  we  wish  to  measure.  It  will  turn  out  that  many  different  aspects  of  this  mea¬ 
surement  converge  at  a  level  of  1  in  104  and,  therefore,  I  shall  not  aim  at  any  higher  precision  in 
the  final  formulas. 

A  second  source  of  complications  is  the  rotational  motions  of  the  reflector.  Any  such  rota¬ 
tion  will  have  the  effect  of  broadening  the  Sun’s  spectral  features  owing  to  the  opposite  Doppler 
shifts  of  the  leading  and  trailing  edges  of  the  satellite.  Of  course,  the  line  spectrum  comes  to  us 
already  broadened  owing  to  the  Sun’s  25-day  rotation,  atmospheric  micro-  and  macro-turbulence, 
other  large  scale  convective  motions,  and  so  on.  Theoretically,  the  additional  broadening  owing 
to  a  satellite’s  rotation  could  be  used  to  compute  the  responsible  angular  velocity.  This  topic  will 
not  be  discussed  herein  except  to  note  that  all  of  the  solar  lines  will  be  equally  broadened  by 
such  motions.  Their  relative  locations  will  not  shift  nor  will  their  absolute  locations.  Hence,  the 
measurement  of  the  satellite’s  center  of  mass  Doppler  shift  is  not  directly  affected  by  a  rotation 
(or  rotations)  about  the  center  of  mass. 

There  are  two  other  reflector-induced  complications.  The  first  has  to  do  with  the  nonuniform 
spectral  reflectivity  of  the  satellite’s  surface  and  the  second  has  to  do  with  its  nonuniform  topog¬ 
raphy.  Because  the  satellite  will  not  be  a  diffuse,  white,  reflecting  sphere,  or  even  a  diffuse,  grey, 
reflecting  solid,  different  portions  of  the  solar  spectrum  will  be  emphasized  after  reflection.  That 
is,  the  gradient  of  the  reflected  continuum  will  not  match  that  of  the  Sun  itself.  Were  we  attempt¬ 
ing  to  utilize  intensity  information,  then  this  would  be  a  severe  problem.  We  are  not,  however, 
and  such  varying  spectral  response  in  the  reflection  coefficient  cannot  alter  the  detected  Doppler 
shift.  Furthermore,  for  practical  reasons  of  instrumental  design,  the  whole  visible  wavelength 
band  from  4000-7000  A  will  not  be  utilized.  Instead,  we  shall  probably  choose  a  several-hundred 
Angstrom  stretch  (or  several  such  disjoint  stretches).  This  will  also  serve  to  minimize  the  effects 
of  a  slowly  varying  background  level  or  a  change  in  the  gradient  of  the  background. 

The  final  reflector-related  complication  comes  about  because  of  topographical  irregularities 
on  the  satellite’s  surface.  These,  coupled  with  the  small  apparent  angular  size  of  the  solar  disc, 
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can  lead  to  specular  reflections.  Such  reflections  are  very  bright  and  of  a  very  short  temporal 
duration.  If  the  radial  velocity  spectrometer  is  using  a  time  averaged  measurement  technique 
(e.g.,  scanning  the  spectrum),  then  the  presence  of  a  very  short  burst  of  very  intense  signal  will 
cause  the  implicit  weighting  of  the  exposure’s  duration  to  be  highly  nonuniform.  This  will  result 
in  a  systematic  error  as  the  epoch  of  the  Doppler  shift  will  not  be  appropriate  owing  to  the  inten¬ 
sity  averaging.  Proper  instrumental  design  and  data  reduction  will  minimize  this  effect. 

Finally,  it  is  necessary  to  discuss  sensor-related  complications.  If  the  spectrometer  is  onboard 
a  spacecraft,  then  the  perceived  Doppler  shift  obviously  includes  a  component  owing  to  the  sen¬ 
sor’s  motion.  As  discussed  above,  this  type  of  complication  is  routinely  handled  in  the  astronomi¬ 
cal  case.  What  makes  the  problem  more  complicated  herein  is  the  use  to  which  the  data  will  be 
put  and  one’s  desire  to  have  the  information  content  simply  relate  to  some  aspect  of  the  dynami¬ 
cal  state  of  the  reflecting  satellite.  If  the  instrument  is  ground-based,  then  this  complication  is 
present  too.  Of  more  importance  may  be  the  additional  corruption  of  the  solar  spectrum  owing 
to  winds  in  the  atmosphere  and  small,  localized  refractive  elements  of  air  that  eventually  build  up 
the  seeing  disc. 

Complications  owing  to  the  Earth’s  atmosphere  are  considerable  but  do  not  affect,  in  princi¬ 
ple,  the  measurement  itself.  Similarly,  a  time  dependent  interplanetary  plasma  will  affect  the  prop¬ 
agation  of  the  light  ray  from  the  Sun  to  the  satellite  and  thence  to  the  sensor.  The  incorporation 
of  these  departures  from  empty  space  and  geodesic  paths  represent  minor  data  reduction  prob¬ 
lems  that  I  acknowledge  but  will  treat  in  depth  after  the  instrument  is  designed. 
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III.  THE  STRUCTURE  OF  THE  ANALYSIS 


It  is  much  more  reliable,  and  esthetically  pleasing,  to  conduct  an  analysis  such  as  this  with 
the  most  complete  physics  ah  initio  and  at  the  most  rigorous  mathematical  level  possible.  Be¬ 
cause  photons  are  inherently  relativistic  objects,  and  gravitational  fields  are  involved,  1  might 
declare  general  relativity  to  be  the  appropriate  venue.  A  little  thought  shows  that  this  represents 
considerable  overkill,  as  the  largest  general  relativistic  effect  is  to  cause  an  additional  redshift 
owing  to  the  gravitational  potential  difference  between  the  target  satellite  or  missile  and  the  sen¬ 
sor.  If  AU  is  the  amplitude  of  the  potential  difference,  then  the  additional  redshift  is  approxi¬ 
mately  equal  to  AU/c.  This  is  —  GM^/cRg  or  of  order  IO"4  km/s.  Hence,  it  can  be  safely  neg¬ 
lected  compared  to  the  desired  precision. 

As  special  relativity  is  the  next  step  down  from  general  relativity,  and  also  represents  the  nat¬ 
ural  covariance  group  of  Maxwell’s  equations,  this  is  the  physical  framework  in  which  the  analy¬ 
sis  shall  be  carried  out.  The  Earth’s  orbital  speed  is  —  30  km/s  or  a  part  in  IO4  (c  =  3  X  105  km/s). 
Similarly,  satellite  or  sensor  speeds  can  not  be  any  higher  (or  much  lower)  in  order  of  magnitude. 
Therefore,  an  eventual  simplification  of  the  fully  covariant  results  to  first  order  will  be  permissi¬ 
ble  and  confluent  with  the  best  attainable  precision  (owing  to  imperfect  knowledge  of  the  fine 
scale  temporal  structure  of  the  Sun,  not  to  the  technique;  it  is  also  true  that  if  we  tried  to 
achieve  centimeters/ second  or  better  precision  that  instrumental  problems  would  rapidly  domi¬ 
nate  the  practical  development).  Finally,  the  last  small  quantity  of  the  problem  is  the  solar  paral¬ 
lax  (i.e.,  the  sensor  to  satellite  vs  sensor  or  satellite  to  solar  distance  ratio).  This  is  —  8 ''8,  or 
again  a  part  in  ten  thousand. 

The  nature  of  the  problem  is  almost  clear  then.  We  commence  with  a  photon  of  wave¬ 
length  X,  frequency  v  emitted  from  the  solar  photosphere.  This  photon  travels  on  a  straight  line 
at  constant  speed  c  until  it  reaches  the  surface  of  some  artificial  satellite  or  ballistic  missile  of 
interest  (the  “target”).  The  photon  is  reflected  from  the  satellite’s  surface  and  travels  toward  the 
sensor.  As  perceived  by  the  sensor,  it  has  frequency  v '  and  wavelength  X'.  By  comparing  X'  to  the 
wavelength  of  a  similar  photon  which  traveled  directly  to  the  sensor,  we  obtain  the  Doppler  shift. 
I  shall  call  the  wavelength  of  this  second  photon  X"  (frequency  v").  There  are  at  least  eight  (8!) 
different  coordinate  systems  hidden  in  the  above  description.  Moreover,  it  is  not  apparent  that  X' 
should  be  compared  with  X"  as  opposed  to,  for  instance,  the  value  of  X  a  sensor  at  the  Earth- 
Moon  barycenter  would  detect,  averaged  over  the  course  of  a  tropical  year.  Once  the  proposed 
analysis  is  fully  worked  out,  modifications  of  this  nature  will  be  relatively  simple  to  tack  on. 


5 


IV.  LORENTZ  TRANSFORMATIONS 


Let  7 r  be  any  contravariant  four-vector.  Let  p0  be  its  time-like  part  and  p^  be  its  space-like 
part,  7r  =  (p0,  p).  An  observer  in  the  inertial  reference  frame  1  measures  n.  Consider  a  different 
observer  in  the  inertial  reference  frame  V.  The  observer  in  1  describes  the  motion  of  the  primed 
reference  frame  as  being  on  a  straight  line  with  constant  velocity  u.  What  will  the  observer  in  V 
detect  for  7r?  The  answer  is  given  by  a  Lorentz  transformation  with  parameter  /?,  £  =  u/c, 

(/?  =|  /?| ),  namely  the  observer  at  rest  in  I'  would  measure  7r'  =  (p'0,  pO  where 

P'ozT(P0-  P/0) 

P]|  =  7(P||  -  /?Po) 

Pi  =  P±  •  (la) 

Herein,  7  “(1  /?2)’,/2  and  the  parallel  and  perpendicular  subscripts  refer  to  components  along 

and  orthogonal  to  u.  As  p  •  u/u  is  the  component  of  p  along  u,  it  follows  that  p  =  p  •  u/u  or 
P||  =  (P  *  u)u/u2.  Consequently,  one  may  write  pi  as  p  -  p  j  and  a  vector  form  of  the  Lorentz 


transformation  as 

7t'=^7(P0  P  /?),  p  +  [(7- UP*  Wu2]u-7p0/?J 

(lb) 

Because  (p  •  u/u2)  u  =  (p  •  /3//J2)  /?,  a  third  form  is 

tt'  =  (t(Po  -  P  '  0)i_P  +  [(7  -  1)  p  \P/P2]P  -  ypod) 

(lc) 

Equations  (1)  and  the  definition  of  an  inertial  frame  are  all  the  special  relativity  that  I  require. 


To  illustrate  Equations  (1),  I  will  consider  the  effects  of  a  Lorentz  transformation  on  two 
different  four-vectors.  One  is  the  space-time  location  vector  p  =  (r0,  r)  =  (ct,  r),  where  t  is  the 
proper  time  andx  is  the  usual  three-dimensional  location,  _r_=  (x,  y,  z).  The  other  is  the  four¬ 
dimensional  wave  vector  k  =  (k0,  k),  where  k0-|_k|  -co/c,  a>  is  the  angular  frequency  of  a  photon 
of  wavenumber  k,  wavelength  A  =  27r/k,  propagating  in  the  direction  of  k_ (frequency  v  =  <d/2tt). 
Note  that  the  four-dimensional  scalar  product  of  p  and  *,  which  is  equal  to  r0k0-_r*  Jc  = 
cot  -v_- _k,  represents  an  invariant,  namely  the  phase  of  the  electromagnetic  disturbance  associated 
with  the  photon. 

If  an  observer  at  rest  in  I  measures  p  and  one  at  rest  relative  to  I'  measures  p',  then 
t'  =  y(t  -n  PIc) 

X  =_L+  [(7  -  ^  yet/?  (2) 

Subsumed  in  Equations  (2)  are  the  two  additional  assumptions  that  at  t  =  t'  =  0  the  origins  of  1 
and  I'  coincided  and  that  no  rotations  are  necessary  to  connect  I  to  I'.  The  arbitrary  translations 
or  rotations  that  might  be  necessary  present  only  mathematical  complexities,  not  any  new  phys¬ 
ics.  From  Equations  (2)  one  can  derive  the  special  relativistic  time-dilation,  the  Fitzgerald-Lorentz 
contraction,  and  so  on. 
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The  derivation  of  the  velocity  transformation  requires  a  bit  more  care  as  v_=  dr/dt  is  not  a 
four-vector.  One  can  add  a  time-like  component  to  _v  and  form  a  proper  four-vector  (this  new 
vector  is  frequently  termed  the  “world  velocity”),  but  I  can  deduce  v'  =  dr'/  dt'  by  a  straight¬ 
forward  differentiation  of  Equations  (2).  The  result  is 

v  +  [(7  -  1)  v  •  u/u21  u  -  yu 

,  ~  .  (3) 

y(l-v^u/c2) 


This  formula  can  be  used  to  compute  relative  velocities.  Consider  two  particles,  labeled  1 
and  2,  with  velocities  vf  and  v2  as  measured  by  an  observer  at  rest  in  some  inertial  frame.  What 
is  the  relative  velocity  of  particle  1  with  respect  to  particle  2?  It  is  the  result  of  applying  Equa¬ 
tion  (3)  to  Vj  with  u  equal  to  v2.  Namely, 


y,  .  —I  +  [(T?.  I)xrv;/V^v2  -Y2^ 

_l  72(>  -Ii  ’  v2/c2) 


(4a) 


where  v'j  is  the  sought-for  relative  velocity,  y2  =  (1  —  y3^)-1 2>  and  J?2  -  v2/c.  As  a  final  point,  a  nice 
form  for  |  v j  |  is 


[I  xi  -  Y21 2  - 1  yi  x  &1 llf¥. 

(1  -v,  •  v2/c2) 


(4b) 


Now  let  us  consider  the  application  of  a  Lorentz  transformation  to  the  four-dimensional 
wave  vector  k  -  (k0,  k).  The  result  is  given  in  Equation  (la), 


k0  =  7(ko  -  k  •  £)>  k||  =  7(k||  -  )3k0),  ki  =  Jcjl 

But  kg  =  cu'/c  and  if  I  write  k  as  kn,  where  n  is  a  unit  vector  in  the  direction  of  propagation, 
then  the  time-like  part  of  k’  =  (kjj,  k^  implies  that 

c o'  -  yaj(l  -  n  •  /}) 

This  embodies  the  Doppler  shift  in  that  co't ^  co  (or  X'?6  X). 


In  order  to  see  the  consequences  of  the  space-like  part  of  the  transformation,  let  6  be  the 
angle  between  u^  and_n.  Then,  k  |  =  k  cos  8  and  kj^  =  k  sin  6.  Writing  kj|  =  k'  cos  6'  and 
k'^  =  k'  sin  6',  it  follows  that 

kj|  =  ky(cos  8  -  (3)  and  kj_  =  k  sin  8 
or,  after  division, 


tan  6 ' 


sin  6 

y(cos  8  - 1 3) 


i.e.,  8¥^8' 


This  shows  that  the  two  observers  perceive  the  light  ray  to  be  propagated  in  different  directions; 
i.e.,  this  formula  embodies  the  aberration  of  light.  Note  that  there  is  aberration  even  if  the  direc¬ 
tion  of  propagation  is  transverse  (0  =  tt/2)  to  the  relative  velocity  vector.  There  is  a  Doppler  shift 
in  this  circumstance  too. 
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V.  FROM  THE  SUN  TO  THE  SENSOR 


Start  a  heliocentric  clock  at  time  t  =  0.  At  this  instant,  the  Sun  has  some  velocity  vector  rel¬ 
ative  to  the  solar  system  barycentric  inertial  reference  frame.  At  t  =  0  construct  an  instantane¬ 
ously  inertial  reference  frame  at  the  Sun’s  center  by  a  translation  of  the  appropriate  amount  and 
via  the  appropriate  Lorentz  transformation.  Also,  at  t  =  0  the  Sun  emits  photons  of  wave- 
number  k,  angular  frequency  co,  isotropically  into  space.  I  now  focus  on  the  photon  which  will  be 
reflected  by  the  target. 

This  reflection  occurs  at  some  time  tT  ( >  0)  later,  as  measured  in  the  original  instantane¬ 
ously  inertial,  heliocentric  reference  frame,  when  the  target  has  a  velocity  vT  as  measured  in  this 
same  reference  frame.  What  wave  vector  does  the  target  perceive?  It  sees  an  incident  four-vector 
*inc,  where  *inc  is  related  to  the  emitted  k  by  a  Lorentz  transformation  with  velocity  _vT.  Thus, 
with  yT  =  (1  -  Vy/c, 

“inc  =  (y-Kko-k-Jj),  k  +  [(7T-  l)k-^r//32r]^T-yTMr)  • 

It  is  understood  that  vT  refers  only  to  the  target’s  velocity  at  time  t  =  tT  as  measured  in  the 
instantaneously  inertial,  heliocentric  reference  frame  of  time  t  =  0. 

Now  consider,  in  the  instantaneously  inertial  reference  frame  of  the  target,  what  the  reflec¬ 
tion  process  looks  like.  If  *inc  has  spatial  part  kinc  and  the  reflected  wave  four-vector  ktc has 
spatial  part  Jcrefl,  then  these  two  three-vectors  lie  in  a  plane  containing  the  outward  surface  nor¬ 
mal  at  the  point  of  reflection;  see  Figure  1.  Also,  the  angle  of  incidence  equals  the  angle  of  re¬ 
flection  in  this  reference  frame.  Thus,  I  shall  be  able  to  compute  as  follows:  I  write  _kjnc  as 


kjnc  “  ^11  +  ^  ± 
krtfi =  -  *ti  *  *J- 


Figure  1.  Reflection  process  in  the  co-moving  reference  frame  of  the  reflector.  The  angle  of  incidence  is 
denoted  by  if/  and  is  also  the  angle  of  reflection  The  outward  surface  normal  is  denoted  by  n  and  the  incoming 
spatial  part  of  the  wave  vector  is  denoted  bykinc.  The  wave  vector  is  decomposed  into  parts  perpendicular  and 
parallel  to  n.  /cy  is  reversed  upon  reflection  so  krej 7  =  -k\\  +  k±. 
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-bine  JSjL  +  JS|| 

where  parallel  and  perpendicular  now  refer  to  the  outward  unit  normal  ji.  Since 
h\\  =  (bine  •  D)  n  and  k±  =  kinc  -  ky 

and  what  happens  upon  reflection  is  that  the  component  orthogonal  to  the  surface  is  reversed 
while  the  component  parallel  to  the  surface  (i.e.,  perpendicular  to  n)  is  unchanged,  all  with  no 
change  of  frequency,  it  follows  that 


ISrefl  —  JL  iS||  l^inc  ^(-inc  - )  — 

Whence, 

*refl  ”  -  —  /^t)’  — inc  “  ^(-inc  *  I?)  JL  j 

(5a) 

with 

kinc  “  k  +  l(yj  “  *)  k  t/Pj]Rt  ~  7Tk0^T 

(5b) 

The  reflected  photon  travels  away  from  the  target.  It  impinges  upon  the  sensor  at  some  time 
ts  (>ty)  as  measured  in  the  original,  instantaneously  inertial,  heliocentric  reference  frame.  At  that 
instant,  in  this  same  reference  frame,  the  sensor  has  velocity  v§.  To  see  what  the  sensor  sees,  I 
must  transform  #creyj,  which  is  expressed  relative  to  the  instantaneously  inertial  reference  frame  of 
the  target  in  Equations  (5),  into  the  sensor’s  reference  frame.  This  implies  performing  a  Lorentz 
transformation  for  velocity  y^'s,  where  v's  is  the  velocity  of  the  sensor  at  t  =  ts  relative  to  what 
the  target’s  velocity  was  at  time  tT.  Equation  (4a)  can  be  used  with  V[  =  vs,  12  =  Iri 


v's  = 


Is  +  [(7t  -  1 )  vs  •  vT/  v^]  vj  -  7XVj 
7-rO  -Is  ■  It/c2) 


P's  =  Is/C 


Applying  Equation  (lc)  to  «ren  I  obtain  the  wave  vector  measured  by  the  sensor, 

«'  =  (7's  (7T<ko  "  Ji-  £t)  ’  tbinc  "  2Qiinc  '  J»)  }  , 

^inc  "  2(kinc  •  n)  n  +  (y's  -  \)  §'s  •  [kinc  -  2(k,nc  •  n)n]0s/(0s)2 
-  TsTt^o  Pj)Psj 


(6) 


(7) 


where  75  =  [1  -  (j8§)2]*w.  The  temporal  part  of  k'  is  oi'/c ,  where  of  is  the  angular  frequency  at 
which  the  sensor  detects  the  photon  originally  emitted  from  the  Sun  with  angular  frequency  co. 
The  Doppler-shift  analysis  involves  examining  a//cy  or  some  similar  quantity.  I  shall  make  the 
reason  for  my  deliberate  fuzziness  clear  below.  Equation  (7)  is  the  penultimate  result. 
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VI.  INTERPRETATION 


Equation  (7)  looks  sufficiently  formidable  that  it  might  be  wise  to  try  and  see  what  it  implies 
in  a  simple  case.  A  two-dimensional  illustration  should  suffice  and  I  take  k  =  k(cos  0,  sin  0,  0). 
Suppose  that  the  reflecting  satellite  is  moving  along  the  x  axis,  =  (/?y,  0,  0),  and  that  the  sen¬ 
sor  is  at  rest,  -0.  Further  suppose  that  the  outwardly  directed  normal  is  in  the  direction  of 
the  negative  x  axis,  _n  =  (-1,  0,  0).  The  other  quantities  that  1  need  are  J3'$  =  -/?T,  y$  =  yT,  and 

7S=  1. 

Note  that  ■  _k  =  kpy  cos  0,  so  Kmc  has  the  form 

K\nc  -  ^7yk(l  _  fir  cos  ®)»  k(7j  -  1)j^t  cos  &T  ~  kyjpj  j 
The  spatial  part  of  this  is  k_inc, 

JSinc  =  k(Tr[C0S  e  -  0t]>  sin  6 *  0) 
and,  as  a  consequence,  1  may  write  Kmc  as 

Kinc  =  k(7r(*  -  h  cos  0)>  7T<cos  0  -  0t),  sin  0,  0  ) 

1  can  also  compute  krefi  =  kinc  -  2(kjnc  •  n)  n, 

k refl  “  k(7rWT  cos  0],  sin  0,  0) 

Whence, 

Krefl  =  k(yj(\  -  (3j  cos  0),  y^pj  -  cos  0),  sin  0,  0  ) 

Putting  all  of  this  into  Equation  (7)  yields  k\ 

k'  =  ky j  ^  1  2/3 j  cos  0  +  (3j,  2(3 j  -  ( 1  +  (3  j)  cos  0,  y  j2  sin  0,  0  ^ 

This  1  need  to  compare  with  k  =  k(l,  cos  0,  sin  0,  0)  =  (k0,  k). 

The  Doppler  shift  in  frequency  is  c(k0  -  k'0)  or 

c(k0  k'0)  =  c  [k  ky  y(l  -  2/?T  cos  0  +  /?y)]  =  2ck/?Tyy(cos  0  -  (3j) 
i.e., 


a>  - 


cu(  1  2(3j  cos  0  +  /?2  ) 


1  -P\ 


If  the  photon  is  normally  incident  (0  =  0)  then  this  takes  a  particularly  simple  form; 


a /(0  —  0)  —  (ju 


—  cu(  1  —  2/3'p)  as  /3y  *  0 


This  should  be  a  familiar  formula  from  monostatic  radar  work.  Astronomers  are  used  to  seeing 
only  one  factor  of  /3T  when  dealing  with  celestial  sources.  That  is  owing  to  the  fact  that  in  such 
work  only  a  single  Lorentz  transformation  is  involved,  i.e.,  into  the  receiver’s  reference  frame.  In 
this  instance  I  have  gone  back  and  forth  yielding  two  factors  of  beta. 
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The  antithesis  of  normal  incidence  is  grazing  incidence  by  a  ray  parallel  to  the  satellite’s  sur¬ 
face.  Now  0  =  7r /  2  and 

co'(0  =  7r/2)  =  a,[(l+/?2T)/(l-/?2T)] 

an  inherently  relativistic  effect  because  it  is  second  order. 


To  see  the  effects  of  aberration,  consider  the  reflection  process  in  the  original  instantaneously 
inertial,  heliocentric  reference  frame.  Then,  Figure  1  looks  like  Figure  2.  While  6  is  both  the 
angle  of  incidence  and  the  angle  between  k  and  (3j ,  O'  is  the  angle  of  reflection  while  n  -  O'  is  the 
angle  between  k'  and  (3j.  Therefore, 


tan(7r  -  O') 


Y-f~  sin  0 

-  [-2/3t  +  ( 1  +  Pj)  cos  0] 


or 


tan  O’ 


(1  -ffftsin  0 
( 1  +  (ij)  cos  0  -  2(3  j 


This  can  be  rewritten  as 


tan  (0'/2)  =  [(1+/?t)/(1  -j8T)]  tan  (0/2) 


I 

t 


Figure  2.  Reflection  process  in  the  inertial  reference  frame  of  the  emitter. 
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Consider  again  the  case  of  normal  incidence;  if  0  =  0,  then  0'  -  0.  This  is  a  simple  statement 
of  the  fact  that  the  wave  vectors  are  always  perpendicular  to  the  wavefronts.  In  the  case  of  graz¬ 
ing  incidence, 

tan  0'  =  -  (1  -  Pj)l (2/3t) 

which  shows  that  for  this  to  occur,  it  must  be  that  /?T<0  when  [0,  tt/2]. 

Now  let  the  sensor  move  too,  /?s  =  (/3S,  0,  0).  The  computation  of  *inc  and  ktc^  is 
unchanged.  From  Equation  (3),  I  deduce  /?'s, 

Whence,  if  krefl  is  the  space-like  part  of  /crefl  [i.e.,  krefl  =  Jc;nc  -  2(kjnc  •  n)  n],  then  the  sensor 
detects  an  original  photon  of  wave  vector  k  as  one  of  wave  vector  k  after  reflection, 

K'  ~  “  Pj  cos  ®)  ~P's  *  iSrefl}’ 

JSrefl  +  [<7s  -  DJs  '  krefl/(/3's)2]^S  "  7'S7T*(  1  "  ^7  cos  9Ws)  • 

A  photon  coming  directly  toward  the  sensor  has  a  perceived  wave  vector  of  k "  [see  Equation  (8)], 
=  (ysk(l  -  /3S  cos  0),  k  +  [(7s  -  1)JS  •  k/02] |j?s  -  7skjSs) 

The  Doppler  shift  c(k'o  -  k'0)  is  given  by 

c(k'o  k'0)  =  ck{7S(l  /3S  cos  0)  -  757^1  - /3T  cos  0  - /3's(^t  _  cos  } 

As  c  -  30  the  Doppler  shift  approaches,  including  all  (32  terms,  2ck(/?T  -  (3S)  (cos  0  -  (3T) 
or 

[1  -  2(3  j  cos  0  +  -  (3c(2(3j  -  cos  0  -  /3s/2)] 

1  -(3\ 

Similarly,  the  aberration  is  given  by 
(7  s)'1  7T1  sin  0 

tan  0'  = - — - - 

(1  -  P’sPj)  cos  0  -  (3  j  +  /3's 

where  (3 =  (/3S  -/3T)/(1  -/3S/3T).  As  c  —  through  all  ( 3 2  terms,  this  reduces  to 

(1  -  (3\ )  sin  0  -  (3J(3q/2  -  /3T)  sin  0 

tan  0'  —  — - t-1 -  s  5 - ^ - 

( 1  +  (3  j)  cos  0  -  2/?t  +  /?s(  1  -  /?T  cos  0) 
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VII.  THE  “DOPPLER”  SHIFT 


A.  THE  REFERENCE  FRAME 

The  analysis  presented  in  Section  V  is  an  instantaneous  one  in  that  it  deals  with  a  particular 
photon  emitted  by  the  Sun  at  a  particular  instant  of  time.  One  might  argue  that  to  close  the 
analysis  the  sensor  should  compare  k  in  Equation  (7)  with  /c",  the  wave  vector  of  a  photon  of 
original  wave  vector  k  emitted  by  the  Sun.  And,  in  Subsection  B  below,  this  is  what  I  shall  do. 
However,  the  reader  should  be  aware  of  the  facts  that  neither  k  nor  k "  are  compelling  choices. 
Indeed,  rather  than  starting  with  k ,  that  is  an  instantaneous  value,  it  might  be  preferable  to  use 
an  averaged  emitted  solar  spectrum.  One  might  average  over  a  month  to  smooth  out  the  lunar 
perturbations  on  the  Earth’s  orbit  about  the  Earth-Moon  barycenter.  Alternatively,  a  longer  time 
span,  such  as  a  year,  might  be  utilized  to  smooth  out  the  reflection  of  the  solar  motion  about  the 
Earth-Moon  barycenter.  Or,  we  might  average  over  a  Jovian  sidereal  period  of  revolution  and 
thereby  smooth  out  the  remaining  major  planetary  perturbations  contributing  to  the  Sun’s  motion 
about  the  solar  system  barycenter.  (I  am  not  neglecting  the  contributions  of  Venus  and  Mars; 
Jupiter’s  period  is  sufficiently  long  that  averaging  over  it  effectively  smooths  their  perturbations 
too.)  No  high  precision  solar  spectrum  has  been  observed  over  a  long-enough  time  interval  to 
provide  the  data  necessary  for  such  averaging. 

Now  consider  the  comparison  for  k.  Should  it  be  k  itself?  If  I  use  the  wave  vector  per¬ 
ceived  by  the  sensor  directly,  then  how  do  we  compute  The  difficulty  is  that,  in  general,  a 
photon  that  left  the  Sun  at  t  =  0,  traveled  along  a  geodesic  to  the  target  satellite,  and  thence 
along  another  geodesic  to  the  sensor,  required  a  different  flight  time  than  a  photon  that  traveled 
directly  to  the  sensor  along  a  third  geodesic.  The  Sun  undergoes  accelerated  motion  so  that 
photons  emitted  at  different  times  with  the  same  wave  vector  relative  to  their  instantaneously 
inertial  heliocentric  reference  frames  have  different  wave  vectors  in  each  others’  reference  frames. 

I  must  either  show  that  the  difference  is  of  the  second  order  for  this  analysis  (as  it  turns  out  to 
be),  correct  for  the  effect,  or  retreat  to  one  of  the  averaged  solar  spectrums  referred  to  above. 

There  are  more  complications.  While  k"  is  the  natural  wave  vector  for  the  sensor  to  use  as  a 
comparison,  assuming  that  a  consistent  reference  frame  for  the  two  ks  can  be  found,  the  sensor 
is  moving  relative  to  the  center  of  mass  of  the  Earth.  The  target  satellite’s  or  missile’s  motion  is 
most  easily  analyzed  in  a  geocentric  reference  frame.  Thus,  a  comparison  of  k "  to  k  does  not 
immediately  yield  a  geocentric  kinematical  variable  (or  combination  thereof).  If  I  throw  the  luni- 
solar  perturbations  of  the  Earth’s  motion  back  into  the  discussion,  then  the  potential  for  compli¬ 
cation  will  increase  further. 

My  principal  purpose  in  recounting  the  above  list  is  to  focus  attention  on  the  nature  of  the 
problem.  The  ultimate  best  coordinate  system  is  not  yet  clear.  The  analysis  contained  herein 
though  is  both  complete  and  rigorous  so  that  any  minor  modifications  can  be  applied  in  a 
straightforward  manner. 
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B.  THE  EXACT  RESULT 


Sloughing  over  for  the  moment  the  exact  description  of  the  reference  frame  of  emission, 
consider  a  photon  with  wave  vector  k  relative  to  the  Sun  that  travels  directly  to  the  sensor. 
Further  suppose  that  it  arrives  at  the  sensor  at  time  ts  as  measured  in  the  original  instantane¬ 
ously  inertial,  heliocentric  reference  frame.  How  will  an  observer  at  rest  with  respect  to  the  sen¬ 
sor  perceive  it?  Because  the  velocity  of  the  sensor,  relative  to  the  original,  instantaneously  inertial, 
heliocentric  reference  frame  is  Vg  at  time  t  =  ts,  the  perceived  wave  vector  k"  is  just  [see  Equa¬ 
tion  (lc)] 

=  (k£,k")  =(ys(k0  -  k  • •  §s),  k  +  [(ys  -  1)  jL‘£s/02s]Js  “  Tsko£s)  <8) 

where  /?s  =  vs/c,  ys  =  (1  -  /Jg)'^. 

C.  THE  POST-NEWTONIAN  APPROXIMATION 

Because  j8  —  1 0~4  for  all  objects  in  all  the  reference  frames  under  consideration,  it  makes 
sense  to  consider  the  limits  of  Equations  (7,8)  as  c  —  In  this  instance,  all  the  gammas 
approach  unity  so 

~  ^  ~  Ji  •  i?T  -  ClSinc  -  2(kinc  ‘  H)ij]  -/S  >  JSinc  "  2(kinc  -  J»)  (k0  "  l &t)£s  ) 

ko-k  £s,  k-k^s) 

and 

k"0  -  k'o  -  -  (fis  -J3j)  ■  k  +  [kinC  -  2(kinc  •  n)  n]  • 

From  Equation  (6)  we  can  see  that 

£s~(k  -  Pj 

which  is  the  Newtonian  result.  From  Equation  (5b)  the  form  for  kinc  reduces  to 
JSinc  ~ ’  ji  _ 

whence 

ko  ”  k'o  —  -2[(^s  -  jij)  •  n]  [k  •  ji]  (9) 

This  result  is  not  transparent,  and  1  now  turn  to  attempting  to  make  it  at  least  translucent. 

D.  THE  CLASSICAL  FORMULA* 

During  the  time  interval  from  t  =  0  until  t  =  tT,  the  incident  photon  travels  from  where  the 
Sun  was  at  t  =  0  in  the  solar  system  barycentric  inertial  reference  frame  to  where  the  target  will 


*  This  form  of  the  analysis  was  suggested  by  E.J.  Kelly  and  was  developed  further  by  R.C.  Raup. 
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be  at  t  -  tT.  The  target’s  geocentric  location  at  time  ty  is  rj,  while  the  Sun’s  geocentric  location 
(a  different  geocenter!)  at  0  is_S0.  Thus, 

ctT=ISo-ly|  =|S(0)  r(tT)| 

Similarly,  the  time  interval  ts  -  ty  represents  the  reflected  photon’s  flight  time  from  the  target  to 
the  sensor.  The  sensor’s  geocentric  (a  third  geocenter)  location  when  the  reflected  photon  arrives 
will  be  ps.  Therefore, 

c(ts  -  tT)  =  \ps  tj\  -  |  p(ts)  -  r(tT)| 

These  two  expressions  can  be  generalized  by  writing  them  as 

ctT  -\$({  ~  ls)  ~L[t  ~  (ts  ~  tj)]l  t=ts  (10a) 

c(ts  -  tT)  =l_p(t)  J1  [t  —  (ts  —  tT)]|  t=ts  •  (I Ob) 

Because  the  sensor,  target,  Earth,  and  Sun  are  all  in  motion,  both  ts  and  ty  are  implicitly 
dependent  upon  solar  system  barycentric  inertial  proper  time  too.  The  Doppler  shifts  induced  by 
these  motions  may  be  obtained  by  differentiating  Equations  (10a)  and  (10b)  with  respect  to  t  and 
then  adding.  Keeping  the  implicit  dependence  of  ts  and  ty  on  t  in  mind,  the  results  are,  after 
evaluation  at  t  =  ts, 

ctT  =  { (1  -  ts)jS0  [1  (ts  '{J?0  ~Ir}/(ctT)  (Ha) 

c(ts  ty)  ={  _ps  -  [I  -  (ts  -  t  j)]  }  *{_PS  “”J[T  }/Cc(tS  “  *t)1  •  (lib) 

With  three  approximations  this  will  reduce  to  the  naive  Newtonian  viewpoint  that  cts  should 
represent  the  instantaneous  time  rate  of  change  of  the  Sun-target-sensor  distance. 

First  note  that  cts  and  cty  are  small,  on  the  order  of  1  in  104  (e.g.,  ~  30  km/s/c).  Hence, 
their  presence  on  the  right-hand  sides  of  Equations  (11)  may  be  ignored.  Second,  observe  that  the 
principal  motion  of  the  Sun  about  the  solar  system  barycenter  is  owing  to  Jupiter’s  attraction. 

This  orbit  has  an  amplitude  of  at  most  M^/Mq  —  5  X  10'3  A.U.  and,  in  fact  —  5  X  10"6  A.U. 

The  period  is  —  12  yr  with  a  linear  tangential  speed  of  —  1  m/s.  Therefore,  since  ty  —  500  s,  1 
can  ignore  the  difference  among _Sq>  Sy,  and _Ss  as  well  as  among  the  corresponding  velocities. 
Third,  note  that  ts  -  tT  —  0.02  s,  so  1  may  approximate  ps  and  ps  by  pT  and  pT.  As  p  and  p  are 
known,  1  can  repair  this  approximation  as  an  additional  term  in  the  planetary  aberration.  Thus, 
the  sum  of  Equations  (11)  can  be  approximately  written,  with  all  now  superfluous  subscripts 
suppressed,  as 

ct  =  <s  -_r)  ■  -|sT7T +  (£  ‘-r) '  ^~rT  '  <l2) 

It  is  also  true  that  p/S  and  r/S  are  —  5  X  10'5.  So,  with  ^  =  S/S  and  only  terms  in  1/S  kept, 
the  final  expression  is  (good  to  1  part  in  104) 

Vj^  =  ct  =  R  +_s-  (S.  l)  +_s_'  (S_— _r)  (L'  s/S)  — _r  *  (S_ -_r) /  S 
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where  R  =  |  p  -  r|  is  the  target-sensor  distance.  Through  the  use  of  s  •  s  =  1  and  the  vector  identity 
(a  X  b)  •  (c  X  d)  =  (a  •  c)  (b  •  d)  -  (a  •  d)  (b  •  c),  this  can  be  recast  as 

Vr  =  R  +  S  •  (S  — i)  +  [(S -i)X  s]  •  (s_X  r/S)  .  (13) 

What  is  the  relationship  between  Equations  (12)  and  (9)?  Equation  (9)  is  the  nonrelativistic 
limit  of  the  fully  covariant  expression  for  kj)  -  kQ.  It  includes  the  unknown  outward  surface  nor¬ 
mal  of  the  target  at  the  time  of  reflection  (e.g.,  ji  at  time  t  =  ty).  It  also  utilizes  velocity  vectors 
referred  to  the  original  (i.e.,  t  =  0),  instantaneously  inertial,  heliocentric  reference  frame.  Sim¬ 
ilarly,  the  space-like  part  of  the  original  wave  vector  also  enters  in  Equation  (9).  However,  in 
Equation  (12)  we  see  geocentric  locations  for  the  target  (r),  the  sensor  (p),  and  the  Sun  (S).  As 
my  earlier  Doppler-shift  results  were  all  expressed  in  the  form  to'  =  oi(I  +  correction),  I  shall 
repeat  this  form  herein. 

The  original  photon  traveled  along  dy,  the  location  of  the  target  at  time  t  =  ty  in  the  instan¬ 
taneously  inertial,  heliocentric  reference  frame  of  t  =  0.  Therefore, 

k  —  k  dy  /  dj 

The  incoming  wave  vector  *jnc,  see  Equation  (5b),  is  a  Lorentz-transformed  version  of  k  =  (k0,  k). 

It  is  turned  into  Krefl,  see  Equation  (5a),  and  then  travels  to  the  sensor.  Thus,  _kref|  travels  from 

d-p  to  ds,  the  location  of  the  sensor  in  the  instantaneously  inertial,  t  =  0  heliocentric  reference 
frame.  It  follows  that  I  can  write  kre<1  as 

krefl  =  krefl(f!s  ~  dy)/|ds  -dy| 

When  I  have  done  this,  to  lowest  order  in  /3,  I  can  derive  an  expression  that  eliminates  n. 

To  start  with, 

K  =  (kQ,  k)  =  k(l,  dT/dT) 

Kinc  r.  k0  -  dj  ■  _@jl dj-,  dTl dT  -j3y) 

Now,  it  follows  that 

Krefl  “*  k  (l  _  lb  ' _fijl dji  dj/ dT  -JTp  -  2[(dT/ dT  -^y)  ‘  j}]  £  j 
=  (k(l  -  dT  -^y/dy),  kref)(dS  -  dy)/|  dS  -  dyl  j 

To  lowest  order, 

krefl  =  k(  1  _^T  ‘  Mj/dj) 

SO 


(1  -  dy  -j8y/dy) 


(js  ~  jy) 

I  ds  -  dy  l 


dy/ dy  -_^y  -  2  [(dy/ dy  -  ^y)  •  Jl]  n^ 


18 


From  Equation  (9)  the  quantity  needed  is  [(/3S  -  0T)  *  n]  (Jc  •  _n)  or  k(n_*  dT/dT)  [(0S  -  (3 j)  *_n].  By 
forming  the  scalar  product  of  the  above  expression  with  /3S  -  (3j  and  keeping  only  the  first-order 
terms,  the  equality 

(E  *  djl dj)  m  “  §j)  *  Jl]  =  Ws  ~  §j)  '  [It/ “  (^s  “  d j)/|  ds  -  djl  ] 

follows. 


The  next  task  is  to  express^  and  d_T  in  terms  of  the  geocentric  quantities^  and_rT.  When 
done  properly,  the  kind  of  implicit  time  dependence  seen  earlier  in  this  Subsection  appears.  As 
my  intent  is  only  to  derive  Equation  (9),  and  all  the  refinements  are  gone  therein,  I  shall  make 
the  necessary  approximations  ab  initio.  Thus, 

d^-p-S^  ,  d  j  =  r  -  Ji 

and 

k'o  =  k  -  k  •  £s  -  2[(J?S  -  lT)  •  n]  (k  •  n) 
so 

(r-S)  (r-S)  R  R. 

c  |_r^— _Sj  c  R 

or 


a)’  =  co(  1  -  V'R/c) 


(14) 


VIII.  CELESTIAL  MECHANICS  IMPLICATIONS 


From  Equation  (14)  it  is  clear  that  the  kinematical  quantity  that  a  radial  velocity  spectrometer 
will  determine  is  Vr.  Looking  at  Equation  (13)  we  see  that  Vr  involves  _r_  and  £_  in  a  complicated 
fashion.  As  £  can  be  written  as_r  =  it,  where  £_  is  a  unit  vector  in  the  direction  of  _r,  it  follows 
that_r  =  r&+  r L  Thus,  since  both  _r_  and  r_  are  involved,  the  determination  of  Vr  provides  infor¬ 
mation  about  the  distance  r,  the  geocentric  radial  velocity  r,  the  geocentric  position  H ,  and  the 

t  ^  " 

geocentric  angular  velocity^.  Clearly,  the  incorporation  of  this  type  of  information  in  a  differen¬ 
tial  correction  procedure  represents  no  new  conceptual  problem  or  difficulty.  On  the  other  hand, 
the  utilization  of  V^in  an  initial  orbit  determination  scheme  presents  unique  difficulties.  A  special 
algorithm  is  required  to  maximize  the  information  content  in  Vr  (which  will  presumably  be  ac¬ 
quired  along  with  the  topocentric  position  and  angular  velocity).  This  algorithm  awaits  invention. 
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